Operating an atom interferometer beyond its linear range 
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Abstract 

In this paper, we show that an atom interferometer inertial sensor, when associated to the 
auxihary measurement of external vibrations, can be operated beyond its hnear range and still 
keep a high acceleration sensitivity. We propose and compare two measurement procedures (fringe 
fitting and nonlinear lock) that can be used to extract the mean phase of the interferometer when 
the interferometer phase fluctuations exceed 2tt. Despite operating in the urban environment of 
inner Paris without any vibration isolation, the use of a low noise seismometer for the measurement 
of ground vibrations allows our atom gravimeter to reach at night a sensitivity as good as 5.5 x 10~'^g 
at 1 s. Robustness of the measurement to large vibration noise is also demonstrated by the ability 
of our gravimeter to operate during an earthquake with excellent sensitivity. Our high repetition 
rate allows for recovering the true low frequency seismic vibrations, ensuring proper averaging. 
Such techniques open new perspectives for applications in other fields, such as navigation and 
geophysics. 
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I. INTRODUCTION 



Atom interferometers [l| are used to develop highly sensitive inertial sensors, which com- 
pete with state of the art "classical" instruments 2l. Applications of such interferometers 
cove. _ field. .0. fundamental pfiy.c. to navigation and .eopl..c. 

For instance, transportable devices are being developed with foreseen applications in the 
fields of navigation, gravity field mapping, detection of underground structures ... 

In most of these experiments, atomic waves are separated and recombined using two- 
photon transitions, induced by a pair of counterpropagating lasers. The inertial force is then 
derived from the measurement of the relative displacement of free-falling atoms with respect 
to the lasers equiphase, which provide a precise ruler. As the inertial phase shift scales 
quadratically with the interrogation time, very high sensitivies can be reached using cold 
atoms along parabolic trajectories [8|,l9|], provided that the experiments are carefully shielded 
from ground vibrations. In the usual geometry where the laser beams are retroreflected on a 
mirror, the position of this mirror sets the position of the lasers equiphase, so that only this 
"reference" optical element is to be shielded from ground vibrations. Such an isolation can be 
realized either with an active stabilization scheme, using a long period superspring 

0,0" 



ul, 



or by using a passive isolation platform 12j. For instance, the use of a superspring allowed 
increasing the interaction time up to 800 ms and reaching a best short term sensitivity to 
acceleration of 8 x 10~^m.s~^ at 1 s An alternative technique, which we study in this 
article, doesn't require any vibration isolation, but exploits an independent measurement of 
ground vibrations, realized by a low noise accelerometer. A technique based on the same 
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principle has already been used with a "classical" corner cube gravimeter 
allowed improving its sensitivity by a factor 7 

In this article, we investigate the limits to the sensitivity of an atomic gravimeter when 
operating without vibration isolation. This transportable gravimeter is developed within the 
frame of the watt balance project led by the Laboratoire National de Metrologie et d'Essais 



(LNE) [16|, Il7l |. We first briefly describe our experimental setup, and recall the usual proce- 



dures for measuring the mean phase of the interferometer. We then introduce and compare 
two measurement schemes (finge fitting and nonlinear lock) that allow operating the sensor 
in the presence of large vibration noise, and show how phase measurements can be performed 
even though the interferometer phase noise amplitude exceeds 27r. These schemes, which 
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use an independent measurement of vibration noise with a low noise seismometer, allow to 
reach good sensitivities without vibration isolation. In particular, we reach a sensitivity as 
good as 5.5 x 10^'^g at 1 s during night measurements, in the urban environment of inner 
Paris. Finally, the robustness of these measurement schemes versus changes in the vibration 
noise is illustrated by the capability of our instrument to operate and measure large ground 
accelerations induced by an earthquake. 



II. LIMITS DUE TO VIBRATION NOISE IN A CONVENTIONAL SETUP 



A. Experimental setup 

The experimental setup, which we briefly recall here, has been described in detail in 
[l^ . About 10'' ^^Rb atoms are first loaded in a 3D-M0T within 50 ms, and further 
cooled down to 2.5 fiK before being dropped in free fall. Before creating the interferometer, a 
narrow vertical velocity distribution of width about 1 cm/s is selected in the \F = 1, mp = 0) 
state, using several microwave and optical Raman pulses. 

The interferometer is then created using Raman transitions js] between the two hyperfine 
levels F = 1 and F = 2 of the ^5*1/2 ground state, which are induced by two vertical and 
counterpropagating laser beams of frequencies uji,uj2 and wavevectors ki,k2. A sequence 
of three Raman pulses (7r/2 — tt — 7r/2) allows to split, redirect and recombine the atomic 
wave packets. The relationship between external and internal state [IJ] allows to measure 
the interferometer phase shift out of a fluorescence measurement of the populations of each 
of the two states. At the output of the interferometer, the transition probability P from 
one hyperfine state to the other is given by P = a + 6cosA$, where 2b is the interferometer 
contrast, and A$, the difference ^ of the atomic phases accumulated along the two paths, 
is given by A$ = —kcs ■ gT"^ [19]. Here k^s = ki — k2 is the effective wave vector (with 
l^effl = ^1+^2 for counter-propagating beams), T is the time interval between two consecutive 
pulses and g is the gravity acceleration. 

The Raman light sources are two extended cavity diode lasers based on the design of 
[2^, which are amplified by two independent tapered amplifiers. Their frequency difference, 
which is phase locked onto a low phase noise microwave reference source, is swept according 
to {uj2 — ^^i){t) = {uj2 — t^i)(0) + at in order to compensate for the gravity induced Doppler 
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shift. This adds aT^ to the interferometer phase shift, which eventually cancels it for a 



B. Conventional measurement procedures 

Maximal sensitivity to phase fluctuations is achieved when operating the interferometer 
at mid fringe, which corresponds to A$ = ±7r/2 . In this case though, variations in the 
offset a can be interpreted as fluctuations of the interferometer phase. A standard technique 
[lo| consists then in recording a full fringe, by measuring the transition probability as a 
function of a controlled phase shift induced on the interferometer. Fitting this fringe then 
allows measuring g. This technique degrades the short term sensitivity as measurements 
performed at the top or bottom of the fringes are not sensitive to phase fluctuations. An 
alternative way consists in using a method inspired by microwave atomic clocks. The phase is 
modulated by ±7r/2 so that the measurement is always performed at mid fringe, alternatively 
to the right and to the left side of the central fringe. From two consecutive measurements 
Pi and Pi+i, the phase error can be estimated. In practice, a correction G x (Pi — Pj+i) 
is added at each cycle to a, in order to stir the chirp rate onto the central fringe. This 
realizes an integrator, whose time constant can be set to a few cycles by adjusting the gain 
G. This locking technique has the advantage of rejecting offset and contrast fluctuations, 
while preserving maximal sensitivity to phase fluctuations. 

C. Influence of vibration noise 

In the case where the duration of the Raman pulses can be neglected, the phase shift A$ 
induced by vibrations is given by 



Jt 

where Zg and Vg are the position and velocity of the experimental setup, and Qs is the 
sensitivity function [21|], given by 



perfect Doppler compensation, for which ao = kes ■ 9- 




(1) 
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The expected sensitivity of the interferometer to g fluctuations of the interferometer is then 



given by a weighted sum of the vibration noise at the harmonics of the cycling rate fc 



211: 
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where crg(r) is the Allan standard deviation of acceleration fluctuations for an averaging 
time r, Sa is the power spectral density of acceleration fluctuations. 

Figure [U displays the power spectral densities of vibrations, measured with a low noise 
seismometer (Guralp CMG-40T, response option 30s) on the platform which is either floating 
(ON) (day time), or put down (OFF) (day time and night time). In the case where the 
platform is OFF, the spectrum is similar to the spectrum measured directly on the ground. 
For our typical parameters, 2T = 100 ms and fc = 3.8 Hz, we calculate using eqJUsensitivities 
at r = 1 s of 2.9 x lO^^g during the day and 1.4 x lO^^g during the night with the platform 
OFF. With the platform ON, the sensitivity is expected to be 7.6 x lO^^g. 
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FIG. 1: Amplitude spectral densities of vibration noise. The black (resp. grey) thick curve displays 
the vibration noise with the isolation platform down (OFF) at day time (resp. night time), while 
the dotted curve displays the vibration noise with the floating platform (ON) at day time. 
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III. VIBRATION NOISE CORRECTION 



A. Correlation between atomic and seismometer signals 

The signal of the seismometer can be used to determine the phase shift of the inter- 
ferometer due to residual vibrations, as measured by the seismometer, 0^^^, which is given 
by: 

{t)vs{t)dt = KffKs / gs{t)Us{t)dt (4) 

JT JT 

where Us is the seismometer voltage (velocity) output and = 4:00.2V/ {m.s~^) is the 
velocity output sensitivity of the seismometer. 

Figure [2] displays the measured transition probability as a function of in the two 
cases of platform ON and OFF, for an interferometer time 2T = 100 ms. The noise is low 
enough in the ON case (fig[2]D)) for the interferometer to operate close to mid fringe, while in 
the OFF case (fig [2^)), interferometer phase noise is larger than 27r, and the interferometer 
signal jumps from one fringe to another. Figure [2] shows the good correlation between 
measured and calculated phase shifts. In the ON case, we find a correlation factor as high 
as 0.94. 




FIG. 2: Correlation between the transition probability of the interferometer and the phase shift 
calculated from the seismometer data, for 2T = 100 ms a) The isolation platform is OFF. Grey 
points : without digital filter, black points : with digital filter, b) The isolation platform is ON. 
Black points : with digital filter. Line : fit to the data, with correlation factor of 0.94. 
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The calculated 0^^^ can thus be used to improve significantly the sensitivity of the mea- 
surement, by applying a post-correction on the transition probability measured at mid fringe. 
This correlation is not perfect though due to the response function of the seismometer, which 
is not flat, and behaves as a low pass fllter with a cut-off frequency of 50 Hz. This response 
function thus limits the efficiency of the vibration rejection. Figure [3] diplays as a continuous 
black line the rejection efficiency as a function of frequency, which is calculated out of the 
seismometer transfer function. 




Frequency (Hz) 



FIG. 3: Efficiency of the vibration rejection as a function of frequency without any processing 
(black straight line), with digital filter (dashed line), with a compensation of a delay of 4.6 ms 
(dotted line). 

B. Digital filtering 

We implemented a numerical filtering of the seismometer signal to compensate for the 
phase lag of the seismometer at intermediate frequencies. The design of the filter is described 
in detail in [12|. It consists in the product of a recursive IIR (Infinite Impulse Response) 
filter, with corner frequencies /o and /i, and a non-causal low-pass filter. The IIR filter 
compensates the phase shift of the seismometer signal and the non causal filter prevents the 
IIR filter from amplifying the intrinsic noise of the seismometer at high frequencies, without 
affecting the phase advance needed to improve the rejection. The total transfer function of 
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the filter is given by 

F(f) = l±lllA I (5) 

where /o, /i and fc are then optimized in order to reach the best sensitivity. This digital 
filtering improves significantly the rejection efficiency, as can be seen in figure [3], where 
it is displayed as a dashed line, for the frequencies /o = 30 Hz, /i = 180 Hz and fc = 
29 Hz. Despite this increase in the rejection efficiency, the gain on the sensitivity, when 
implementing this filter in the ON mode, was limited to 25% only [l^, which we attributed 
to excess noise of the seismometer arising from coupling between horizontal and vertical 
axes. 



C. Cross couplings 



In order to detect these couplings, we recorded simultaneously the seismometer outputs 
along the three directions, calculated three corrections, one along each axis (only the vertical 
correction was numerically filtered though) and fitted the transition probability measured at 
mid fringe with a linear combination of the three corrections. The result of this fit showed 
couplings of 4% and 5% with the horizontal axes. We finally determined the influence of 
these couplings onto the sensitivity of the measurement, by comparing the Allan standard 
deviation of the phase fluctuations in the case where the correction is performed only with 
the vertical correction (ID), or with the optimal combination of the three (3D). The results 
are shown on figure HI where the sensitivity is expressed relatively to g. Using the three 
corrections allows to remove a bump that appears when using the ID correction. This 
indicates that horizontal vibration noise, as it also appears in the vertical seismometer 
signal, adds noise when performing a ID correction. 



D. Efficiency of the filter without vibration isolation 

The digital filter is much more efficient in the OFF mode, as one can see on figure [2] a) 
where the noise on the interferometer fringes is significantly reduced when seismometer data 
are processed with the digital filter. In that case, the dominant contribution of the vibration 
noise onto the degradation of the sensitivity corresponds to frequencies around 10 Hz, for 
which the effect of the filter improves the rejection efficiency from 10 dB to about 30 dB 
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FIG. 4: Sensitivity to g with ID and 3D corrections. The measurement was realized during the 
day, with a floating platform. 

12|. 

E. Case of a pure delay 

We later noticed that the phase lag of the seismometer signal varies almost linearly with 
respect to frequency in the 1-100 Hz band, with a slope corresponding to a delay of about 
5 ms. The phase shift of the seismometer can thus be compensated for, by simply shifting 
the acquisition of the seismometer data by this delay. We measured the correlation factor 
as a function of the delay, with the platform OFF, and found an optimal delay of 4.6 ms. 
The rejection efficiency for this optimal delay is displayed as a dotted line on figure [3l 
Surprisingly, we find a correlation similar to the optimal digital filter, despite a significantly 
different behavior of the rejection efficiency versus frequency. 
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IV. MEASUREMENTS PROTOCOLS 



A. Standard procedures 

The standard measurement protocols described above need that phase fluctuations remain 
significantly smaller than 2tt. This requires to reduce the interferometer duration in the 
OFF mode to 2T < 20 ms. For 2T = 20 ms, the integrator scheme described above 
allows reaching sensitivities of 1 x lO^^g at 1 s when applying no correction to the measured 
transition probability, of 5 x lO^^g when correcting without filtering, and of 1.5 x lO^^g when 
correcting with digital filtering. The simple post-correction (without filter) thus improves 
the sensitivity by a factor 2, and the digital filter improves it further by a factor 3.5. Better 
performances are expected with large interrogation time for which the transfer function 
of the interferometer filters more efficiently high frequency vibration noise. In order to 
operate the interferometer with large interrogation times despite excess noise, we propose 
two alternative measurement procedures described in the following subsections. Both are 
based on the combination of measurements of the transition probability and of (p^-j^ by the 
seismometer. Though developed for the case of large vibration noise, these techniques can 
be extended to low vibration noise by adding a well controlled phase modulation. 

B. Fringe fitting 

The first technique simply consists in fitting fringes, as in |l^, except that here the phase 
of the interferometer is now scanned randomly by vibration noise. The signal displayed in 
figure [2] and obtained when plotting the transition probability versus 0^^^, calculated with 
the digital filter, can be fitted by the function P = a + 6cos(r70fj^ + 6(f)), where a,b,ri and 
6(j) are free parameters. Due to the infiuence of the seismometer transfer function, rj will in 
general differ from 1. In practice, we operate the interferometer close to the central fringe, 
which corresponds to a small phase error 6(f). Every 20 points, we perform a fit of the 
signal and extract a value for the phase error 6(j)m- We then calculate the Allan standard 
deviation of the 6(j)m in order to determine the sensitivity of the measurement. Note that 
this fitting procedure is not very efficient if the noise amplitude is significantly less than 
2tt, because the interferometer signal remains close to the bottom of the central fringe. An 
additional and perfectly controlled phase modulation of ±7r/2 is thus apphed in order to 
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optimize the sensitivity of the interferometer to phase fluctuations. Moreover, the sensitivity 
improves by about 50% when taking cross couphngs of the seismometer into account, which 
can be reahzed by adjusting the data with a hnear combination of the corrections along 
three directions T^rjjCpl-^ -, where j = x,y,z and 4>vib,j the phase shift calculated out of the 
filtered seismometer data along axis j. 

C. Nonlinear lock 

The lock procedure described in IIIBI can be adapted in the case where the phase noise 
exceeds 27r. Let's consider the measurement at cycle i of the transition probability Pj 

Pi = a — b cos{{kcsg — a)T^ + Si) = a — 6(cos e cos Si — sin e sin Si) (6) 

where e = {Ksg — a)T'^ is the phase error and Si is the phase shift induced by residual 
vibrations, estimated from the seismometer signal. We assume here that the phase error e 
varies slowly, so that we can consider it as constant between three consecutive measurements. 
Eliminating a and cos e from the following three equations 

Pi-i = a — 6(cos Si^i cos e — sin S'j_i sin e) 

Pi = a — 6(cos Si cos e — sin Si sin e) 
Pi+i = a — 6(cos Si+i cos e — sin Si+i sin e) 

gives 

bBi sin e = Ai 

with 

Ai = (cos^i+i - cosS'i)(i^_i - Pi) - {cosSi-i - cos Si){Pi+i - Pi) 

Bi = (cos 5*4+1 — cosS'j)(sinS'j_i — sin 6*4) — (cosS'j_i — cos S'j)(sin S'j+i — sins'!) 

In order to stir the chirp rate onto the doppler shift rate, an iterative correction is applied 
on a according to 

ai+2 = ttj+i + K '^^'2 A (7) 
1 + i^i 

where is a positive gain. Here is used as a pseudo inverse of bBi with b ^ 1/2, in 
order to prevent the correction from diverging when Bi is close to zero. Choosing K < 
guarantees the stability of the servo loop. 
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D. Adaptation of the nonlinear lock 



When phase fluctuations are signiflcantly less than 1 radian, Bi becomes much smaUer 
than 1 (note that Bi is null in the absence of vibration noise, which implies that the lock 
scheme doesn't work, as it is not able to stir the chirp rate), so that is not a good pseudo- 
inverse of Bi. This decreases the effective gain of the loop, which can be compensated for 
either by increasing K, or by replacing j^hi with ^txr^-, where cr^ is the standard deviation 



The scheme is then modifled by adding extra phase shifts in order to increase the sen- 
sitivity to phase fluctuations. A simple phase modulation of ±7r/2, which implies that the 
interferometer operates alternatively at the right and left sides of the central fringe, is not 
sufficient, as in that case, Bi is still null for null vibration noise. With a 3-phases modu- 
lation (— 7r/2, 0, 7r/2), Bi = 1 for null vibration noise, and replacing j+W ^'^^^ +b^+b'^ ' 
with B the mean of BiS, guarantees the full efficiency of the lock, whatever the amplitude 
of vibration noise. 

The lock technique can be further modifled to flrst determine and servo the vibration 
phase coefficients rjj. The phase of the interferometer is e + Si + S(l)i, where is a con- 
trolled additional phase shift (alternatively — 7r/2, 0, 7r/2), and the vibration phase Si is (best 

approximated by) 'Erij(f)'^-j^ j, where j = x,y,z and (p^nij is the phase shift calculated out of the 
seismometer data along axis j. At the i-th measurement. Si is calculated by Y^^=i Vj,i4'vib,j,i^ 
where rjj^i — rjj — Srjj^i. Pi is thus given by 



of the BiS. 
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Pi^a- bcos{S(j)i + J2 Vj,i<l>vib,j,i + e + XI H,i(f>vib,j,i) 



3 




where Si = 5(f)i + E^Li Vj,i<Pvib,j,i 



Generalizing the algebra above, one gets 
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(8) 
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where 



Cj^i = (cos Si+i - cos 5,)(</)^ifej. sin S^-l - 0^^^^. ^ sin Si) 
-(cos Si_i - cos S^) sin S^+i - ^^^^j^^ sin Si) 

Chirp rates and vibration phase coefficients are then corrected according to 
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where Lj is the gain for direction j. Such nonhnear feedback and estimation algorithms are 
inspired from Lyapounov stabihty theory, for the main loop given by eq. Ul and adaptive 



techniques, for the estimation of parameters rjj (see 22| for a tutorial presentation of such 
techniques and I23I] for a more advanced one). 

Figure [5] displays the evolution of the vibration phase coefficients during a two day mea- 
surement. The time constant of the lock is about 200 s (see inset). Note that the vertical 
phase coefficient 773 differs significantly from 1, and is different at day and night times, which 
can be attributed to a change in the vibration noise PSD. Moreover, the lock converges to- 
wards horizontal phase coefficients of about 5%, in agreement with the values previously 
determined with the fit. 



E. Comparison of the two techniques 

Figure [6] displays the Allan standard deviation of g fluctuations for 2T = 100 ms, with 
the two techniques described above (fringe fitting and nonlinear lock), during day and night 
times. The vibration phase shifts were calculated out of the 3D signals, using the optimal 
delay of 4.6 ms. We obtain equivalent sensitivities at 1 s of 2.7 x 10~^g (resp. 1.8 x 10~^g) 
with the nonlinear lock (resp. fringe fitting) technique during the day, and 8.5 x lO^^g 
(resp. 5.5 x lO^^g) during the night. We find that the fit of the fringes is slightly better 
than the lock technique, by about 50%. The efficiency in removing vibration noise from the 
gravimeter signal can be calculated from the ratio of the sensitivities obtained here with 
the calculated contribution of the vibration noise (see section Hi Cp . A gain from 11 to 25 is 
obtained depending on the technique and noise conditions. 
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FIG. 5: Evolution of the vibration phase coefficients, during a measurement reahzed using the 
nonhnear lock scheme, with initial settings rjj^Q = (0,0,1). The graph on the left (resp. right) 
displays the vertical (resp. horizontal) phase coefficient (s). 

Best sensitivities are obtained during night measurements, as the vibration noise in the 
1-10 Hz band is significantly lower. We reach at best an equivalent sensitivity as low as 
5.5 X lO^^g at 1 s when fitting f ring es, which is only 4 times worse than our best reported 



value with the platform floating 12|], and only twice larger than the sensitivity obtained in 
our laboratory with a commercial FG-5 corner cube gravimeter in the same vibration 
noise conditions. 

These two techniques were also compared in a numerical simulation, where the phase of 
the interferometer was generated randomly as the sum of two independent terms = 0i + 02, 
with Gaussian distribution of standard deviations (Xi and o"2. (pi simulates the vibration 
phase noise measured by the seismometer 0^^^, and 02 the phase difference between the real 
vibration phase noise and 0^^^,. We then implemented the two techniques with such simulated 
data, with a2 = 0.02 rad and with ai ranging from 0.06 to 30 rad. For each technique, we find 
the corresponding sensitivity of the interferometer at 1 shot 0"$ and calculate a normalized 
sensitivity by dividing 0"$ with a2- We verified that this normalized sensitivity does not 
depend on a2- The results of the simulations are displayed in figure [7] and for both techniques 
the normalized sensitivity exhibits the same behavior. It increases for vibration noise larger 
than a few hundreds mrad, for which linear approximation of the transition probability is no 
longer valid, and finally saturates for large vibration noise. This degradation is due to the 
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FIG. 6: Allan standard deviation of g fluctuations versus averaging time. Measurements with the 
nonlinear lock technique at day (resp. at night) are displayed as full stars (resp. open stars). 
Measurements with the fringe fitting technique at day (resp. at night) are displayed as full circles 
(resp. open circles). 

non linearity of the transition probability versus interferometer phase: measurements at top 
and bottom of the fringes have no sensitivity to phase fluctuations. The simulation confirms 
that this degradation is higher for the lock technique than for the fringe fitting technique, 
as observed in the measurements. In particular, for (T2 = 3 rad, which corresponds roughly 
to day conditions, we find normalized sensitivities of 1.28 and 1.80 for the fringe fitting and 
lock techniques. The ratio of the sensitivities is thus 1.4, in reasonable agreement with the 
measurements. 

F. Investigation of systematic effects 

It is important to verify that the techniques presented here provide an accurate mea- 
surement of the interferometer phase, free from any bias. The lock procedure, which is 
intrinsically non-linear, could in principle induce such a bias. The numerical simulation 
indicates that none of the two techniques suffer from such systematics. This was confirmed 
experimentally by performing differential measurements, alternating the standard integra- 
tion technique described in III Bl with the lock procedure described in[IVDl in the case where 
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FIG. 7: Numerical simulation of the normalized sensitivity of the interferometer as a function of 
the vibration noise standard deviation. Black squares (resp. open circles) display the sensitivity 
degradation for the nonlinear lock (resp. fringe fitting) technique. 

the platform was ON and thus the noise level low. The difference between the two techniques 
was found to be 0.3 ± 0.8/iGal, which is consistent with no bias. Moreover, the two tech- 
niques were compared together during the day with the platform OFF, which corresponds 
to a noise level of a^s = 3 rad. The difference for a 6 hours measurement was found to be 
—2 ± 4 mrad, which corresponds to —5 ± 10/iGal, which is also consistent with no bias. 



G. Interest of the nonlinear lock procedure 

The main advantage of the nonlinear lock scheme is a better time resolution. Indeed, the 
time constant of the lock loop can be reduced to a few cycles only, so that a time constant 
< 1 s can be reached. In comparison, fitting the fringes requires to fit data by packets of 
at least 20 cycles for optimal sensitivity, which reduces the time resolution to about 5 s. 
Both techniques can operate with low vibration noise. Indeed, the fit of the fringes can 
also be adapted by modifying the phase modulation to add measurements performed at 
the top and bottom of the interferometer, in order to constrain the sinusoidal fit (doing so, 
sensitivity will as well be degraded because these measurements are not sensitive to phase 
fluctuations). We finally illustrate the efficiency of the lock algorithm by demonstrating its 
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robustness versus large changes in the vibration noise. Figure [8] displays the measurement 
during an earthquake of magnitude 7.7 that occurred in China on the 20th of March 2008. 
The gravimeter detects efficiently the occurrence of seismic waves, of period about 20 s. 
As our seismometer, of long period 30 s only, measures these vibrations with a large phase 
lag of about 1 radian, they are not efficiently removed from the gravimeter phase shift by 
the lock algorithm. They thus appear as a clear and well resolved signal in the gravimeter 
data. This demonstrates the robustness of our system versus large excitations, which is 
not the case for traditional absolute corner-cube gravimeters, which have neither adequate 
repetition rate (usually about 0.1 Hz) nor sufficient dynamic range, due to the finite range 
of the superspring mechanism. Note that the use of a longer period seismometer would in 
principle allow removing these low frequency vibrations from the gravimeter data. 




FIG. 8: Fluctuations of the gravimeter signal during the earthquake of magnitude 7.7 that occurred 
in China on March 20*^^ 2008. Data were obtained with the nonlinear lock procedure. 



V. CONCLUSION 

In this paper, we demonstrate that an atom interferometer can reach high sensitivities 
without vibration isolation, when using an independent measurement of vibrations by a low 
noise seismometer. We develop here several measurement protocols that allow determining 
the mean phase of the interferometer, even when the interferometer phase noise amplitude 
exceeds 27r. In particular, fitting the fringes scanned by vibration noise allows reaching a 
sensitivity as low as 5.5 x lO^^g at 1 s during night measurements. This performance is 
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obtained with a rather short interaction time (2T = 100 ms), for which the vertical length 
of the interferometer corresponds to a few centimeters only. 

The techniques presented here are of interest for the realization of a portable atom 
gravimeter, with potential application to geophysics and gravity measurements in noisy 
environments. A compact gravimeter, associated with a good AC accelerometer and oper- 
ating at a high repetition rate would reach fairly high sensitivities, without much hardware 
isolation against ground vibrations. Moreover, in contrast with other classical instruments, 
such as ballistic corner cube gravimeters, a high sensitivity would still be reached in the 
presence of earthquakes, if using a long period seismometer (100 s) to measure vibration 
noise. 

More generally, these techniques can be extended to differential measurements with atom 
interferometers, such as gradiometers and cold atom gyroscopes. In particular, the phase 
difference can easily be extracted from the fits of the two interference patterns. Strong 
interests of these techniques lie in the ability of extending the dynamic range of the sensors, 
and of extracting the inertial phase without bias. 
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